Kramers degeneracy in a magnetic field 
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o 
o 



(N 



Revaz Ramazashvili* 
ENS, LPTMS, UMR8626, Bat. 100, Universite Paris-Sud, 91405 Orsay, France 

(Dated: March 21, 2009) 

In this article, I analyze the symmetries and degeneracies of electron eigenstates in a commen- 
surate collinear antiferromagnet. In a magnetic field transverse to the staggered magnetization, a 
hidden anti-unitary symmetry protects double degeneracy of the Bloch eigenstates at a special set 
of momenta. In addition to this 'Kramers degeneracy' subset, the manifold of momenta, labeling 
the doubly degenerate Bloch states in the Brillouin zone, may also contain an 'accidental degenera- 
cy' subset, that is not protected by symmetry and that may change its shape under perturbation. 
These degeneracies give rise to a substantial momentum dependence of the transverse g-factor in 
the Zeeman coupling, turning the latter into a spin-orbit interaction. 

I discuss a number of materials, where Zeeman spin-orbit coupling is likely to be present, and 
outline the simplest properties and experimental consequences of this interaction, that may be 
relevant to systems from chromium to borocarbides, cuprates, hexaborides, iron pnictides, as well 
as organic and heavy fermion conductors. 

PACS numbers: 75.50.Ee 
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I. INTRODUCTION 

Antiferromagnetism is widespread in materials with 
interesting electron properties. Chromium [1] and its 
alloys [2, 3], numerous borocarbides [4], electron- and 
hole-doped cuprates [5, 6], iron pnictides [7], various or- 
ganic [8] and heavy fermion [9-11] compounds all have an 
antiferromagnetic state present in their phase diagram. 
The physics of these antiferromagnetic phases has been 
a subject of active research. 

In this article, I study the response of electron Bloch 
eigenstates in an antiferromagnet to a weak magnetic 
field. I concentrate on the simplest case: a centrosym- 
metric doubly commensurate collinear antiferromagnet, 
shown schematically in Fig. 1, where the magnetization 
density at any point in space is parallel or anti-parallel 
to a single fixed direction n of the staggered magnetiza- 
tion, and changes sign upon primitive translation of the 
underlying lattice. 

In a paramagnet, the double degeneracy of the Bloch 
eigenstates is commonly attributed to symmetry under 
time reversal 9 - and, indeed, perturbations that break 
time reversal symmetry (such as ferromagnetism or a 
magnetic field) do tend to remove the degeneracy. Yet 
violation of 8 alone does not preclude degeneracy: in a 
commensurate centrosymmetric Neel antiferromagnet, as 
in a paramagnet, all Bloch eigenstates enjoy a Kramers 
degeneracy [12] in spite of time reversal symmetry being 
broken in the former, but not in the latter. 

In an antiferromagnet, the staggered magnetization 
sets a special direction n in electron spin space, ma- 
king it anisotropic. A magnetic field along n removes 
the degeneracy of all Bloch eigenstates, as it does in a 
paramagnet. By contrast, in a transverse field, a hidden 
anti-unitary symmetry protects the Kramers degeneracy 



of Bloch eigenstates at a special set of momenta. 

Generally, in d dimensions, the manifold of momenta, 
corresponding to doubly degenerate Bloch states in a 
transverse field is (d — l)-dimensional; within a subset 
of this manifold, the degeneracy is dictated by symmet- 
ry. This is in marked contrast with what happens in a 
paramagnet, where an arbitrary magnetic field lifts the 
degeneracy of all Bloch eigenstates. For brevity, in this 
article I often refer to the manifold of momenta, labeling 
the degenerate Bloch states in the Brillouin zone, as to 
the 'degeneracy manifold'. 

As a consequence of the Kramers degeneracy of the 
special Bloch states in a transverse field, the transverse 
component g± of the electron g-tensor vanishes for such 
states. Not being identically equal to zero, g± must, 
therefore, carry a substantial momentum dependence, 
and the Zeeman coupling Hzso must take the form 

Hzso = -/is [5|| (H|, • a) + g±(p)(H ± ■ tr)] , (1) 

where H|| = (H • n)n and Hi = H H|| are the longi- 
tudinal and transverse components of the magnetic field 
with respect to the unit vector n of the staggered mag- 
netization, [Ib is the Bohr magneton, while g» and gj_(p) 
are the longitudinal and transverse components of the 
g-tensor. 

This significant momentum dependence of g± (p) turns 
the common Zeeman coupling into a kind of spin-orbit 
interaction Hzso whose appearance and key pro- 
perties are the focus of this work. Zeeman spin-orbit 
coupling may manifest itself spectacularly in a number 
of ways, which will be mentioned below and discussed in 
detail elsewhere. 

The symmetry properties of wave functions in mag- 
netic crystals have been studied by Dimmock and 
Wheeler [13], who pointed out, among other things, that 
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magnetism not only lifts degeneracies by obviously low- 
ering the symmetry but also may introduce new ones. 
This may happen at the magnetic Brillouin zone (MBZ) 
boundary, under the necessary condition that the mag- 
netic unit cell be larger than the paramagnetic one [13]. 

For a Neel antiferromagnct on a lattice of square sym- 
metry, the response of the electron states to a magnetic 
field was studied in Ref. [14] using symmetry arguments, 
and in Ref. [15] within a weak coupling model. The 
present work is a detailed presentation of recent results 
[16]. It revisits Ref. [14], extends it to an arbitrary crys- 
tal symmetry and to a finite as opposed to infinitesimal 
magnetic field, and uncovers a rich interplay between the 
symmetry of magnetic structure and that of the under- 
lying crystal lattice. At the same time, the present work 
extends Ref. [13] by allowing for an external magnetic 
field - to show how, at special momenta, the Kramers 
degeneracy in an antiferromagnet may persist even in a 
transverse magnetic field. 

This work treats antiferromagnetic order as static, ne- 
glecting both its classical and quantum fluctuations. This 
excludes from consideration strongly fluctuating antifer- 
romagnetic states such as those near a continuous phase 
transition, be it a finitc-tcmpcraturc Neel transition or a 
quantum (T = 0) critical point. At the same time, the 
single-electron Bloch eigenstates, considered hereafter, 
must be well-defined. As in a normal Fermi liquid state, 
this does not rule out strong interaction between elec- 
trons, but simply requires temperatures well below the 
Fermi energy. Finally, to justify the neglect of quantum 
fluctuations, the ordered magnetic moment must be of 
the order of or greater than the Bohr magneton. 

As a consequence, the present theory applies to anti- 
ferromagnets (i) deep inside a commensurate long-range 
antiferromagnetic state, and far enough from any con- 
tinuous Neel transition, finite-temperature or quantum, 
(ii) with an ordered moment noticeable on the scale of the 
Bohr magneton, and (Hi) far below both the Neel and the 
effective Fermi temperatures. All materials mentioned in 
Section IV are meant to be considered under these con- 
ditions. 

The article is organized as follows. Section II opens 
with a reminder of how, in spite of broken time re- 
versal symmetry, all Bloch eigenstates in a commensu- 
rate collinear antiferromagnet retain Kramers degenera- 
cy, provided there is an inversion center [12]. Then I 
show how, even in a transverse magnetic field, a hid- 
den symmetry of antiferromagnetic order may protect 
the Kramers degeneracy for certain Bloch states. 

Section III establishes several properties of the single- 
electron spectrum in a weakly coupled antiferromagnet, 
subject to a transverse magnetic field. 

Section IV contains the analysis of simple examples, 
that may be relevant to specific materials from chromium 
to organic conductors, from borocarbides to underdoped 
cuprates and to various heavy fermion metals. The Dis- 



cussion section reviews the findings, and examines them 
in the light of earlier work, while the Appendices present 
various technical details. 



II. GENERAL ARGUMENTS 

It is convenient to begin by describing the crystal host 
symmetry in the absence of magnetism, with the average 
magnetization density notionally set to zero [13]. I re- 
fer to this as to the paramagnetic state symmetry, even 
though the symmetry of the actual paramagnetic state 
may be different, for instance due to a lattice distortion 
upon transition. Unitary symmetries of the paramag- 
netic state form a group, ft, which includes a set of ele- 
mentary translations T a by primitive translation vectors 
a. Time reversal 9 being indeed a symmetry of the para- 
magnetic state, the full symmetry group g of the para- 
magnetic state includes ft, and products of 9 with each 
element of ft: g = h+9-h; put otherwise, ft is an invariant 
unitary subgroup of g. 

Antiferromagnetic order couples to the electron spin 
<x via the exchange term (A r • er), where A r is propor- 
tional to the average microscopic magnetization at point 
r. In keeping with the arguments of the Introduction, 
fluctuations of A r are neglected. Being of relativistic 
origin, spin-orbit couplings of the crystal lattice to the 
electron spin and to the magnetization density are also 
neglected. This is a good approximation in a broad range 
of problems, at the very least at temperatures above the 
scale set by the spin-orbit coupling (see Section V for de- 
tails). In this 'exchange symmetry' approximation [17], 
magnetization density and electron spin are assigned to 
a separate space, independent of the real space of the 
crystal; this makes coordinate rotations and other point 
symmetries inert with respect to A r and er. 

A nonzero A r changes sign under time reversal 9, and 
removes the symmetry under primitive translations T a , 
thus reducing the symmetry with respect to that of para- 
magnetic state. In a doubly commensurate collinear an- 
tiferromagnet, A r changes sign upon T a : A r+a = — A r , 
while T a leaves A r intact: A r+2a = A r . Even though 
neither 9 nor T a remain a symmetry, their product 9T a 
does (see Fig. 1). In a system with inversion center, so 
does 9T a X, where X is inversion. The importance of the 
combined symmetry 9T a X will become clear in the next 
subsection. 

Together with the uniaxial character (A r at any point 
r pointing along or against the single direction n of stag- 
gered magnetization), these relations define a commen- 
surate collinear Neel antiferromagnct via transformation 
properties of its microscopic magnetization density. 
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FIG. 1: (color online). Doubly commensurate collinear anti- 
ferromagnet on a simple rectangular lattice. In the absence 
of magnetism, time reversal and primitive translations T x 
and T y , shown by dashed arrows, are symmetry operations. 
In the antiferromagnetic state, neither of the three remains 
a symmetry, but the products 9T X and 0T y , shown by solid 
arrows, do, as illustrated by filled spin arrows. Small dashed 
rectangle at the center is the Wigner-Seitz cell boundary in 
the paramagnetic state, while the shaded hexagon is its an- 
tiferromagnetic counterpart. Notice that neither of the point 
group operations interchanges the two sublattices, hence any 
point symmetry of the lattice, including inversion T, remains 
a symmetry of the antiferromagnetic state. 



Kramers degeneracy in zero field 

The combined anti-unitary symmetry 6*T a Z gives rise 
to a Kramers degeneracy [12]: If |p) is a Bloch eigenstate 
at momentum p, then #T a Z|p) is degenerate with |p). 
Since 9 and X both invert the momentum, both |p) and 
6>T a Z|p) carry the same momentum label p. Formally, 
this is verified by the action of any translation Tb, that 
remains a symmetry of the antiferromagnetic state: 

T b #T a X|p) = 0T a T b Z|p) = 0T a JT_ b |p) = 
= #T a Ze- ! P- b |p) =e*P- b 0T a Z|p). (2) 

At the same time, |p) and 6*T a Z|p) are orthogonal. This 
follows from Eqn. (22) of Appendix A as soon as one 
chooses O = T a J, = |p), and \<f>) = T a Z6»|p). Recall- 
ing that (T a Z) 2 = -6> 2 = 1, and hence (9T a l) 2 = -1, 
one finds 



<p|JT a %) = -(p|2T a 0|p>. 



(3) 



Thus, in spite of broken time reversal symmetry, in a 
centrosymmetric commensurate Neel antiferromagnet all 
Bloch states retain a Kramers degeneracy. 

Kramers degeneracy in a transverse field 

Generally, a magnetic field H lifts this degeneracy. 
However, in a transverse field, a hidden anti-unitary sym- 



metry may protect the degeneracy at a special set of 
points in the Brillouin zone, as I show below. 

In an antiferromagnet, subject to a magnetic field H, 
the single-electron Hamiltonian takes the form 



H = H + (A r • a) - (H • a) 



(4) 



where the 'paramagnetic' part 7io is invariant under in- 
dependent action of T a and 9, and g^is is set to unity. In 
the absence of the field, all Bloch eigenstates of Hamil- 
tonian (4) are doubly degenerate by virtue of Eqn. (3). 
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1 = n xh 



A r+a =-A 



FIG. 2: (color online). Relative orientation of A r , A r + a , Hy 
and Hi. To see the combined symmetries of Table I, notice 
that 6 flips both A r and H, while T a leaves H intact, but 
inverts A r . Unit vectors n, h and 1 = n x h are defined as 
shown. 

Consider the symmetries of Hamiltonian (4) , involving 
a combination of an elementary translation T a , time re- 
versal 9, or a spin rotation U m (</>) around an axis m by 
an angle <j>. These symmetries are listed in Table I; the 
relative orientation of A r , Hii and Hj_ is shown in Fig. 
2. The transverse field Hj_ breaks the symmetries JJ n (4>) 
and T a #, but preserves their combination at 4> = tt, i.e. 
U n (7r)T a 6'. Acting on an exact Bloch state |p) at mo- 
mentum p, this combined anti-unitary operator creates a 
degenerate partner eigenstate U n (7r)T a 6'|p), which is or- 
thogonal to |p) everywhere in the Brillouin zone, unless 
(p • a) is an integer multiple of ir (in other words, unless 
p lies at a paramagnetic Brillouin zone boundary): 



<p|U n (7r)T a %) = e - 2 *(P' a >(p|U n (^)T a %). 



(5) 



Equation (5) follows from Eqn. (22) of Appendix A for 
O = U n (^)T a , |V) = |p>, and \<f>) = U n (^)T a %) as 
soon as one observes, that [U n (7r)T a 6 1 ] 2 = T a = T-2 a . 
In a magnetic field, double translation T2 a remains a 
symmetry; according to the Bloch theorem, it acts on 
|p) as per T 2a |p) = e 2i ( p ' a )|p), thus leading to (5). 

Notice, however, that the eigenstate U n (7r)T a 6'|p) car- 
ries momentum label — p rather than p. By contrast with 
the case of zero field, combining U n (7r)T a with inver- 
sion X no longer helps to produce a degenerate partner 
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TABLE I: Symmetries of a collinear doubly commensurate 
antiferromagnet. The left column is for zero field, central col- 
umn for a transverse magnetic field Hi , and the right column 
for a longitudinal field Hii. As above, T a denotes elementary 
translation by a. As shown in Fig. 2, unit vector n is collinear 
with A r , unit vector h points along Hi, and unit vector 1 is 
defined via 1 = n x h. For a general orientation of the field 
(H || 7^ 0, Hi 7^ 0), not shown in this table, \J\(tv)9 is the 
only surviving symmetry, where U m (</>) denotes spin rotation 
by angle <f) around the m axis. 



(A r ■ a) 


(A r • a) + (Hi • er) 


(A r • a) + (H,| - a) 


Un(<« 




Un«>) 


U h (7T)T a 


U h (7T)T a 




Ui(7r)T a 






T a 6> 


Un(7T)T a 




U h (7r)0 




U h (7T)0 


Ui(tt)0 
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eigenstate at the original momentum p: since 9, U n (7r), 
and T a T all commute, and since [XU n (7r)T a #] = 1, 
equation (22) of Appendix A for O = XU n (7r)T a , 
tp = |p), and <j) = ZU n (7r)T a (9|p) only confirms, that 
(p|ZU n (7r)T a #|p) equals itself. 

Thus, for an exact Bloch state |p) at momentum p, 
the anti-unitary symmetry U n (7r)T a produces an or- 
thogonal degenerate eigenstate U n (7r)T a 6>|p) at momen- 
tum —p. The two momenta p and — p are different, with 
one key exception. It occurs for p at the magnetic Bril- 
louin zone boundary, given a unitary symmetry U, that 
transforms — p into a momentum, equivalent to p up to a 
reciprocal lattice vector Q of the antiferromagnetic state 
[13]: 



-Up = p + Q. 



(6) 



In this case, the eigenstate WU n (7r)T a #jp) carries mo- 
mentum label p + Q = p, is degenerate with |p) and or- 
thogonal to it, thus explicitly demonstrating Kramers de- 
generacy at momentum p in a transverse field. This re- 
sult is general: combined with any momentum-inverting 
anti-unitary symmetry, equation (6) leads to Kramers 
degeneracy at momentum p [13]. The simplest illustra- 
tion, where U is the unity operator, is given by p = Q/2 
and shown in Figs. 5 and 6(a) for two particular cases. 
These and other examples are described in Section IV. 
Notice that, at p = Q/2 (with U = 1), the degeneracy 
in a transverse field is guaranteed even for a low crystal 
symmetry, provided an inversion center. 

Also notice, once more, that U, X and other point sym- 
metries above are inert with respect to spin as a conse- 
quence of the 'exchange symmetry' approximation [17]. 

Hamiltonian (4) and the subsequent analysis ignored 



A r (Hi) , ; L A (Hi) 




Ap(H ) ( 



1 = n xh 



r+a 



(Hi) 



A r+a (Hi) 



FIG. 3: (color online). Relative orientation of A r (_Hi) and 
A r + a (-Hi), tilted by a transverse field Hi. The component 
A r L (ffi) along the field is periodic under elementary trans- 
lation T a : A^fffi) = A.^r(H±); the component along the 
zero-field A r is anti-periodic: Aj! +a (i/i) = -A r (iix). 



the response of the antiferromagnetic order to the trans- 
verse field Hj_. This, however, does not affect the set of 
points, where Kramers degeneracy in a transverse field 
is protected by the anti-unitary symmetry WU n (7r)T a 6>. 
Upon application of Hj^ , the Neel sublattices tilt towards 
the field, making it convenient to present A r as 



(7) 



where A^(H±) points along Hj_, and A}(H±) points 
along n, as shown in Fig. 3. Since A^r +SL (H±) = A^r(H±) 
and Al +EL (Hj_) = —AI(Hj_), the second column of Table 
I remains intact upon replacing A r in Hamiltonian (4) 
by A r (H x ) ofEqn. (7). 

If the antiferromagnetic unit cell is a multiple of its 
paramagnetic counterpart, the magnetic Brillouin zone 
boundary contains a set of points, that do not belong to 
the paramagnetic Brillouin zone boundary (for example, 
see Figs. 5 and 6). In the paramagnetic state, no two 
points of this set, separated by antiferromagnetic recip- 
rocal lattice vector Q and satisfying condition (6), can be 
declared equivalent. As a curious consequence, the mag- 
netic group of such a wave vector is not a subgroup of its 
paramagnetic counterpart [13]. Hence the degeneracy, if 
present, does hinge on magnetic order. 



III. CLUES FROM WEAK COUPLING 

Additional insight into the locus of states, that remain 
degenerate in a transverse magnetic field, is afforded by 
a weak-coupling single-electron Hamiltonian in a doubly 
commensurate collinear antiferromagnet. Let Q be the 
antiferromagnetic ordering wave vector (see the examples 
below); A r creates a matrix element (A -<r) between the 
Bloch states at momenta p and p + Q (for simplicity, I 
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neglect its possible dependence on p). Sublattice canting 
in a transverse field is taken into account in Appendix B. 
In magnetic field H, and at weak coupling, Hamiltonian 
(4) takes the form [15] 



U 



6 P - (H • a) (A • cr) 



(A • cr) 6 P+Q - (H • <r) 



(8) 



where e p and e P +Q are single-particle energies of TCq in 
(4) at momenta p and p + Q, and the 'bare' g-tensor in 
(H • er) is omitted for brevity. 

In a purely transverse field Hi, this Hamiltonian can 
be diagonalized simply by choosing the z-axis in spin 
space along H^, and the i-axis along A. As a re- 
sult, Hamiltonian (8) splits into two decoupled pieces: 
7Yi(p, HjJ for the amplitudes |p; ]) and |p + Q; J,), and 
7Y 2 (p, Hj_) for the amplitudes |p; J.) and |p + Q; ]): 



Wi ( 2)(p,Hl) = 



A e p+Q ± H± 



(9) 



The spectra £1(2) (p) of Ti-i(2) are given by 



£ 1(2) (p,H ± )=, ?p ± v /A 2 + [C P T(H 



(10) 



with (Hj_ -cr) = H± corresponding to Hi, and (Hj_ -er) — 



—H± to H2, and with r] p = 



_ e p +e p+ Q 



, and C P 



2 ' ""^ — 2 

The same spectrum can be obtained by excluding, say, 
|p + Q;tr) from the eigenvalue equation for (8), but it 
is important to keep in mind that cr in (10) no longer 
describes spin, but rather pscudospin: since (Hj_ -cr) docs 
not commute with the Hamiltonian, the eigenstates of 
T~ti(2) are superpositions of spin-up and spin-down states. 

Equation (10) illustrates a number of points. Firstly, 
the electron spectrum acquires a gap of size 2A. Second- 
ly, in the absence of magnetic field, each eigenstate is 
indeed doubly degenerate, in agreement with the argu- 
ments, encapsulated in Eqn. (3). Thirdly, Eqn. (10) 
shows, that the degeneracy persists in a transverse field 
(and, therefore, ,gj_(p) in Eqn. (1) vanishes) whenever 
C p = 0. Barring a special situation, this equation defines 
a surface in three dimensions, a line in two, and a set of 
points in one. This result for the dimensionality of the 
manifold of degenerate states hinges solely on the sym- 
metry of the antiferromagnetic state and holds beyond 
weak coupling, as shown in Appendix C. Furthermore, 
as shown above, this manifold must contain all points, 
satisfying Eqn. (6): the points, where the degeneracy is 
enforced by symmetry. Finally, expansion of Eqn. (10) 
to first order in (Hj_ • cr) yields the expression for g±(p) 
in (1) within the weak coupling model (8): 



flx(p) 



C P 



A 2 + C^ 



(11) 



At the end of the preceding Subsection, I showed that 
tilting of the Neel sublattices in a transverse field does not 
affect the set of points, where Kramers degeneracy in a 
transverse field is protected by the anti-unitary symmetry 
U\J n (n)T a 9. However, generally, the rest of the degen- 
eracy manifold is not protected by symmetry, and may 
change shape upon crystal deformation, or under ano- 
ther perturbation. For instance, while leaving intact the 
symmetry-protected set of degeneracy points, the sub- 
lattice canting may change the shape of the degeneracy 
manifold 5_l(p) = compared with C P = 0. This effect 
is discussed in Appendix B. 

Put otherwise, the degeneracy manifold may be di- 
vided into two parts. The first part is the 'Kramers de- 
generacy' subset of special momenta, fixed by conspiracy 
between the anti-unitary symmetry U n (7r)T a ^ and the 
crystal symmetry. This 'Kramers' subset is insensitive to 
perturbations that leave intact the crystal symmetry of 
the material. The rest is an 'accidental' degeneracy sub- 
set, whose geometry, by contrast, may vary under pertur- 
bations, that do not affect the crystal symmetry, but only 
alter the microscopic parameters of the system. The di- 
vision of the degeneracy manifold into the 'Kramers' and 
the 'accidental' degeneracy subsets is well illustrated by 
the examples of two-dimensional rectangular and square 
symmetry antiferromagncts in Section IV. 



Spectral symmetries in momentum space 

The spectrum of Hamiltonian (8) enjoys a number 
of symmetries. Firstly, inversion symmetry makes the 
spectrum even under inversion. At the same time, 
5_l(p)(Hj_ • cr) must also be even under inversion, which 
implies 



S-l(-p) = flj.(p)- 



(12) 



The anti-unitary symmetry TJ\(tt)9 in the last line of Ta- 
ble I is another reason for g± (p) to be even under inver- 
sion, as U^Tr)^ turns g_i_(p)(Hj_ -tr) into gj_(— p)(Hj_-<r). 

Periodicity doubling due to antiferromagnetism mani- 
fests itself more interestingly. The ordering wave vec- 
tor Q being a reciprocal lattice vector, any Bloch eigen- 
state at momentum p must have a degenerate part- 
ner eigenstate at momentum p + Q. Usually, this im- 
plies Q-pcriodicity of any given band: e(p) = e(p + Q), 
which is the case, for instance, in a longitudinal field 
H || , where £(p) undergoes the common Zccman split- 
ting £(p) — > £(p) ± Hn. In a Neel antiferromagnct in 
a transverse field, this is not the case: for a general p, 
£i(2)(p + Q) t^£i(2) (P)- Instead, 



£i(p + Q ) H L )=£,(p,H L ) ) 



(13) 



while both £i(p) and £ 2 (p) of Eqn. (10) are invariant 
under momentum shift p — > p + 2Q. These properties 
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are illustrated in Fig. 4, showing the splitting of a one- 
dimensional conduction band in a transverse field. 




FIG. 4: (color online). One-dimensional conduction band 
splitting, (a) Conduction band ['+' sign in front of the square 
root in Eqn. (10)], split by a transverse field. Here Eqn. (13) 
is illustrated by spectrum (10) for a nearest-neighbor hopping 
e p = tcosp in Hamiltonian (8). Notice that, in spite of pe- 
riod doubling in real space, neither of the two split spectra 
has period it; instead, both are 27r-periodic, but map onto 
each other upon translation by ir. Also notice that the two 
split sub-bands intersect at the symmetry-enforced degener- 
acy point p = ±~, as they should (see the one-dimensional 
example of the following section), (b) Conduction band, split 
by a longitudinal field. By contrast with (a), each split sub- 
band is 7r-periodic, and degeneracy is lifted for all momenta. 

The reason behind (13) is that, as long as H± ^ 0, 
neither Hi nor H2 in Eqn. (9) is invariant under the mo- 
mentum boost p — ► p + Q, in spite of the Hamiltonians 
(4) and (8) both having doubled periodicity. Rather, 

«i(p + Q,Hl)=W 2 (p,Hi), (14) 

which is made explicit by subsequent exchange of the 
diagonal matrix elements, and leads to Eqn. (13). 

Zeeman splitting corresponds to the difference 
£i(p, Hj_) — £2(p, Hj_); hence it changes sign upon mo- 
mentum shift by Q. At the same time, direct in- 
spection shows that Hamiltonians Hi and H2 in Eqn. 
(9) turn into one another upon inversion of Hj_: 
Wi(p, — HjJ =Ti-2{p, HjJ. Combining this with Eqn. 
(14), one finds that momentum boost by Q accompanied 
by inversion of Hj^ is a symmetry of both Hi and Ti.2- 

Wi(2)(p + Q,-Hx)=Wi (2 )(p,Hl). (15) 

For the transverse Zeeman term g±(p)(H±-cr), this yields 

g±(p + Q) = -g±(p)- (16) 

Combined, Eqns. (12) and (16) lead to 

<7±(|+P) = -S±(§-P). (17) 

This implies not only that g±{p) must vanish at p = 
but also that g±(% + p) is an odd function of p. The 
conclusions of this subsection hold after the sublattice 
tilting is taken into account (see Appendix B). 



IV. EXAMPLES 

In this section, I describe the manifolds of degenerate 
states for a number of concrete examples and thus show, 
that the Zeeman spin-orbit coupling (1) is at work in 
many materials of great interest. It gives rise to various 
interesting phenomena, some of which are outlined in the 
subsection 'Experimental signatures' of Section V. 

Kramers degeneracy in a transverse field and the Zee- 
man spin-orbit coupling (1) will manifest themselves 
whenever carriers are present at or near the manifold 
of degenerate states gj_(p) = 0. In a weakly doped an- 
tiferromagnetic insulator, this will happen whenever the 
relevant band extremum falls at or near the manifold of 
degenerate states. In an antiferromagnctic metal, this oc- 
curs when the Fermi surface crosses this manifold. Hence, 
for metals, I mention the Fermi surface geometry when- 
ever known. 

Between these two limiting cases of a weakly doped an- 
tiferromagnetic insulator and an antiferromagnetic metal 
with a large Fermi surface, the experimental manifesta- 
tions of the Zeeman spin-orbit coupling will be quanti- 
tatively different. On top of this, certain effects will be 
sensitive to the geometry of the degeneracy manifold and 
its intersection with the Fermi surface, as well as the ori- 
entation of the staggered magnetization with respect to 
the crystal axes. A detailed discussion of these effects 
will be presented elsewhere. 

When selecting the examples below, the preference 
was given to materials, available in high-purity samples, 
where de Haas-van Alphen oscillations were observed, 
and where magnetic structure was unambiguously cha- 
racterized by neutron scattering. As explained in the 
Introduction, the results of this work apply to materials 
well inside a long-range antiferromagnetic phase, and far 
enough from any critical point, quantum or classical. For 
both quantum and thermal fluctuations of antiferromag- 
netic order to be negligible, the ordered moment shall be 
noticeable on the scale of the Bohr magneton, and the 
sample shall be kept well below both the Neel and the 
effective Fermi temperatures. 

One dimension 

In one dimension, the magnetic Brillouin zone bounda- 
ry reduces to two points p = ±§S which in fact coincide 
up to the antiferromagnetic wave vector Q = — , that is 
also the reciprocal lattice vector of the antiferromagnetic 
state (see Fig. 5). In terms of the general condition (6), 
this is the simplest case: U = 1. 

As a result, at p = ±|-, the two exact Bloch states in 

1 1 2a 1 

a transverse field, |p) and 0T a U n (7r)|p), correspond to 
the same momentum p, and are degenerate by virtue of 
6>T a U n (7r) being a symmetry. Equation (5) guarantees 
their orthogonality, thus protecting Kramers degeneracy 



at momentum p = against transverse magnetic field. 
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FIG. 5: (color online). The paramagnetic (p = ±— ), and 
the antiferromagnetic (p = ±5^) Brillouin zone boundaries 
of a one-dimensional Neel antiferromagnet. In the antiferro- 
magnetic state, the two points p = i^- are identical up to 
the reciprocal lattice vector Q = — of the antiferromagnetic 
state. At these two points, anti-unitary symmetry U n (7r)T a 6> 
protects the Kramers degeneracy against transverse magnetic 
field. 



Two dimensions, rectangular and square symmetries 

Now consider a two-dimensional antiferromagnet on a 
lattice of rectangular or square symmetry, with the or- 
dering wave vector Q = (71", n). In a transverse magnetic 
field, the degeneracy persists on a line in the Brillouin 
zone, by virtue of Eqn. (10). I will show that, in the rect- 
angular case, the degeneracy line must contain the point 
£ at the center of the magnetic Brillouin zone (MBZ) 
boundary (i.e. the star of p = Q/2), shown in Fig. 6(a). 
In the square symmetry case, the degeneracy persists at 
the entire MBZ boundary (Fig. 6(b)). The MBZ in Fig. 
6 is the reciprocal space counterpart of the Wigner-Seitz 
cell of the magnetic state (Fig. 1), and the ordering wave 
vector Q = (n, n) connects points X and Y in Figs. 6(a) 
and (b). 

Consider a Bloch state |p) at momentum p in a trans- 
verse field. As discussed in Section II, the eigenstate 
#T a U n (7r)|p) at momentum — p is degenerate with |p) 
and, according to Eqn. (5), must be orthogonal to it 
unless (p • a) is an integer multiple of it - put other- 
wise, unless p belongs to the paramagnetic Brillouin zone 
boundary. At points E, X, and Y, momenta p and — p 
coincide up to a reciprocal lattice vector of the antifer- 
romagnetic state. However, at points X and Y (as well 
as at the entire vertical segment of the MBZ boundary 
in Fig. 6(a)), (p • a) is an integer multiple of n; hence 
|p) and #T a U n (7r)|p) are not obliged to be orthogonal 
there as per Eqn. (5). Thus, £ is the only point at 
the MBZ boundary, where the two degenerate states |p) 
and 6>T a U n (7r)|p) are orthogonal and correspond to the 
same momentum: Dashed arrows in figure 6(a) show, 
that, for a generic point p' at the MBZ boundary, no 
symmetry operation relates — p' to a vector, equivalent 
to p'. Therefore, it is only at point E, that the sym- 




(a) 



(b) 



FIG. 6: (color online). Geometry of the problem, (a) The 
Brillouin zone for a simple rectangular lattice (the rectan- 
gle), and its antiferromagnetic counterpart (MBZ, shaded 
hexagon). Thick (red) curve, passing through point E, shows 
a typical degeneracy line <?x(p) = 0. At the MBZ boundary, 
only momentum p at point E is equivalent to — p up to a 
reciprocal lattice vector of the antiferromagnetic state; for a 
generic p', shown by the dashed arrow, this is not true, (b) 
The Brillouin zone of a simple square lattice and its antiferro- 
magnetic counterpart (shaded diagonal square). The degen- 
eracy line must contain the entire MBZ boundary, shown in 
red. Point T is the Brillouin zone center, points X and Y 
lie at the centers of the paramagnetic Brillouin Zone edges. 
Point E lies at the center of the MBZ boundary. 



metry protects Kramers degeneracy against transverse 
magnetic field. As in the one-dimensional example of 
the previous subsection, in terms of Eqn. (6) this corre- 
sponds to the simplest case of U = 1. 

This can be illustrated by a nearest-neighbor hopping 
spectrum 



e p = t [cos + rj cos p y ] 



(18) 



in the weak-coupling example of the previous subsection: 
for rectangular symmetry (?/ 7^ 1), spectrum (10) in a 
transverse field remains degenerate at a thick (red) line, 
sketched in Fig. 6(a). Upon variation of rj ^ 1, the line 
changes its shape, but remains pinned at the star of wave 
vector p = Q/2 (i.e. at point E) in Fig. 6(a). In terms 
of the preceding subsection, the star of p = Q/2 is the 
'Kramers' subset of the degeneracy manifold, while the 
rest of the degeneracy line in Fig. 6(a) is the 'accidental' 
degeneracy subset. 

Promotion from the rectangular symmetry (?) ^ 1) to 
that of a square (rj = 1) brings along invariance under 
reflections a\ 2 in cither of the two diagonal axes 1 and 
2, passing through point T in Fig. 6(b). As a result, 
the eigenstate eri(9T a U n (7r)|p) at momentum CT2P (Fig. 
6(b)) is also degenerate with |p) and orthogonal to it, 
as one can show analogously to the examples above. In 
terms of the general condition (6), this means hi = <J\.2- 



A momentum p at the MBZ boundary in Fig. 6(b) 
differs from 02P by a reciprocal lattice vector; thus the 
two momenta coincide in the nomenclature of the an- 
tiferromagnetic Brillouin zone. Hence, for a square- 
symmetry lattice in a transverse field, the degeneracy 
is of a 'Kramers' (i.e. symmetry-protected) nature at the 
entire MBZ boundary as shown in Fig. 6(b). In this 
case, barring a particularly pathological band structure, 
the degeneracy manifold is exhausted by its 'Kramers' 
subset. 

In accordance with the symmetry arguments above, 
for the toy nearest-neighbor hopping spectrum (18) at 
the square symmetry point 77 = X, the degeneracy line of 
Eqn. (10) coincides with the MBZ boundary as shown 
in Fig. 6(b). By contrast, for rectangular symmetry, it 
is Eqn. (10) that restricts the degeneracy in a transverse 
field to a line in momentum space, and it is the symmetry 
that pins this line at point £ at the middle of the MBZ 
boundary, as shown in Fig. 6(a). 

Now, g±(p) can be expanded in a vicinity of the de- 
generacy line g±(p) = 0. With the exception of higher- 
symmetry points, such as point X in Fig. 6(b), the lead- 
ing term of the expansion is linear in momentum devia- 
tion Sp from the degeneracy line: 



S-l(p) 



Hp • Sp 



(19) 



where 3 p /7i is the momentum gradient of <?j_(p) at point 
p on the degeneracy line. As mentioned in the previ- 
ous Section, inversion symmetry makes g±(p) even un- 
der inversion. Therefore, H p changes sign upon inversion, 
which is consistent with Eqns. (13) and (16), that require 
<7_i_(p) to change sign upon momentum shift by Q. 

As shown in the last subsection of Section III, g±(p) is 
an odd function of the deviation Sp from point £ (the star 
of p = Q/2) in Figs. 6(a) and (b). Therefore, expansion 
of <?_i_(p) around point £ cannot contain an even power 
of Sp. 
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FIG. 7: (color online). Schematic drawing of the commensu- 
rate magnetic structure of antiferromagnetic chromium. 




FIG. 8: (color online). The Brillouin zone and the Fermi 
surface of commensurate antiferromagnetic chromium, (a) 
Dashed lines show the paramagnetic Brillouin zone bounda- 
ry. The solid cube inside it is the magnetic Brillouin zone 
(MBZ), with high symmetry points indicated. In the anti- 
ferromagnetic state, point H is equivalent to point V at the 
center of the Brillouin zone (not shown), (b) The Fermi sur- 
face sketch of paramagnetic chromium. As in (a), dashed lines 
show the paramagnetic Brillouin zone boundary, and the solid 
cube inside it is the MBZ. The nearly octahedral hole Fermi 
surface is centered at point H, and nearly spherical electron 
'balls' are located at face centers X of the MBZ boundary. To- 
gether with the nearly octahedral electron surface, centered 
at point r (not shown), these electron balls form the electron 
'jack'. A set of hole ellipsoids is centered at points N in the 
middle of the magnetic Brillouin zone edges. 



This subsection is devoted to commensurate antiferro- 
magnetism in chromium - the simplest of magnetic or- 
ders, occurring in this textbook spin density wave metal. 
Chromium crystallizes into a b.c.c. lattice, and under- 
goes various magnetic and structural transitions upon 
variation of temperature, pressure, or alloying [1-3]. 

Below the Neel temperature TV of about 311 K at am- 
bient pressure, chromium develops weakly incommensu- 
rate antiferromagnetism with ordered moment of about 
0.5 hb per atom at 4.2K. However, strain - or doping 
with some 0.1 to 0.3% of a transition metal (such as Mn, 
Re, Rh, Ru, Ir, Os or Pt [2]) - eliminate incommensu- 
rability in favor of commensurate order with wave vec- 
tor [001], shown in Fig. 7. Commensurate order has 



also been observed and much studied in thin films of 
chromium, often with an enhanced Neel temperature and 
ordered moment [18]. This article neglects fluctuations 
of magnetic order and, conveniently, the high Neel tem- 
perature of chromium facilitates experimental access to 
T -C Tjv, where thermal fluctuations are suppressed. 

The paramagnetic and the antiferromagnetic Bril- 
louin zones for bulk commensurate antiferromagnetic 
chromium are shown in Fig. 8(a). An arbitrary mo- 
mentum at the MBZ boundary becomes equivalent to 
its opposite upon reflection in a properly chosen plane. 
Similarly to the two-dimensional square-symmetry exam- 
ple above, this equivalence is up to a primitive wave vec- 
tor of the antiferromagnetic reciprocal space. Hence, in 
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TABLE II: Simple properties of some of the studied samples 
of Celri3, UGa3 and UI113: the Neel temperature Tjv, the 
ordered magnetic moment M, the Sommerfeld coefficient 7, 
the residual resistivity po, and the residual resistivity ratio 
p(300A-)/po. 



Material 


T N 
(K) 


M 
(/xb) 


7 

(mJ/K 2 -mol) 


Po 
(/iS7cm) 


p(300K) 
PO 


Celn 3 


10.1 


0.5/Ce 


130 


0.5 [20] 


35 [20] 






[19] 


[20] 


0.6 [21] 


100 [22] 


UGa 3 


64 


0.75/U 


52 


1.2 [26] 


38 [25] 




[23] 


[23] 


[24-26] 




81 [26] 


UIn 3 


88 


1/U [27] 


50 [24] 


0.66 [28] 


130 [28] 



a transverse magnetic field, the Kramers degeneracy sur- 
vives at the entire magnetic Brillouin zone boundary in 
Fig. 8(a). 

The disappearance of g± (p) affects electrons at two dif- 
ferent sheets of the Fermi surface, sketched in Fig. 8(b): 
those at the nearly spherical electron parts, centered at 
points X in the middle of each MBZ face, and those at 
the hole ellipsoids, centered at points N in the middle 
of each MBZ edge. For the former, the leading term of 
the expansion is linear in the momentum deviation Sp± 
from the flat face of the MBZ boundary. For the lat- 
ter, the leading term of the expansion is quadratic near 
each MBZ edge, since <7_i_(p) vanishes at each of the two 
intersecting faces of the MBZ boundary. 

Celn 3 , Ulna, UGa 3 ... 

A number of cerium and uranium binary intermetallics 
of simple cubic Cu 3 Au structure, such as Celn 3 , CeTl 3 , 
Ufn 3 , UGa 3 , UT1 3 and UPb 3 , turn antiferromagnetic at 
low temperatures. High-purity samples of Celn 3 , UIn 3 , 
and UGa 3 have made it possible to characterize magnetic 
order and electron properties of these materials rather 
comprehensively. Some of the basic properties of the 
samples are shown in Table II. 

At low temperature, all three develop a type-II anti- 
ferromagnetic structure with wave vector Q = [555], 
shown in Fig. 9(a) for Celn 3 . The materials remain 
normal metals down to the lowest temperatures probed, 
with the Sommerfeld coefficient substantially enhanced 
by comparison with that of a simple metal (see the fourth 
column of Table II versus about 0.65 mJ/K 2 -mol for Ag). 

Of the three materials, Celn 3 has been scrutinized the 
most. Its early studies were driven by interest in valence 




FIG. 9: (color online). Geometry of Celn 3 in real and in re- 
ciprocal space, (a) Cubic unit cell of Celn 3 , showing Ce atoms 
and their magnetic moments. Indium atoms (not shown) are 
positioned at the face centers of the unit cell, (b) Cubic Bril- 
louin zone of paramagnetic Celn 3 and, inside, its antiferro- 
magnetic counterpart. Darker shading marks the degeneracy 
surface. 



[29] and magnetic [19] fluctuations, in the nature of its 
magnetic order [19], in large mass enhancement [30] and 
related questions. Subsequent research focused on the re- 
duction of Tat under pressure, and on superconductivity, 
discovered near the critical pressure p c , where the Neel 
temperature is about to vanish - as well as on marked de- 
parture from Landau Fermi liquid behavior, found in the 
normal state near p c [22, 31, 32]. The most recent work 
included de Haas-van Alphen oscillation measurements 
[21, 33], electron-positron annihilation experiments [34], 
and interpretation of the former [35] . 

According to Fig. 9(b), the Magnetic Brillouin Zone of 
the three metals enjoys full cubic symmetry. Its square 
faces belong to the paramagnetic Brillouin zone bound- 
ary (p • a) = ±7r, where Eqn. (5) does not enforce de- 
generacy; however, g±(p) does vanish at the hexagonal 
MBZ faces, marked by darker shading in Fig. 9(b). 

According to de Haas-van Alphen measurements [21, 
36] and to calculations [37], one sheet of the Fermi sur- 
face of Celn 3 is nearly spherical, and has radius of about 
— where a is the lattice constant. Hence this sheet 

a 2 1 

comes close to the point L in Fig. 9(b), which is the very 
same distance — ^ away from the Brillouin zone center. 
Disappearance of g±(p) necessarily affects the dynamics 
of an electron on this sheet in a transverse field. 

Near a generic point at an MBZ face, far from its edges, 
leading terms of the expansion of g±(p) are linear in 
transverse deviation of momentum from the MBZ face 
as per Eqn. (19), with H p normal to the MBZ boundary. 
Near the edges, joining the neighboring hexagonal faces 
in Fig. 9(b) - for instance, near the points E and W - 
the leading terms become quadratic. 
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Uranium nitride 

Uranium nitride (UN) presents another example of in- 
terest. This heavy fermion metal has a face-centered cu- 
bic lattice of NaCl type, shown in Fig. 10(a). Below 53K, 
it develops type-I antiferromagnetic order, with ordered 
moment of about 0.75/is per uranium atom [38], and the 
Sommerfeld coefficient of 50 mJ/K 2 -mol [39]. The Neeel 
temperature of UN drops under pressure, vanishing at 
about 3.5 GPa. Recent experiments [39] studied the low- 
temperature resistivity near the critical pressure on sam- 
ples with residual resistivity po of about 2.3 /ificm, and 
the residual resistivity ratio p(300K) / po of the order of 
10 2 . 




FIG. 10: (color online). Geometry of uranium nitride (UN) 
in real and in reciprocal space, (a) F.c.c. cubic unit cell of 
UN, showing U atoms and their magnetic moment orientation. 
Nitrogen atoms are shown by open circles, (b) Dashed lines 
define the Brillouin zone boundary of paramagnetic UN; the 
square prism inside it is the antiferromagnetic Brillouin zone. 
Its entire boundary defines the degeneracy surface <?_i_(p) = 0. 

The real-space sketch of magnetic structure of UN is 
shown in Fig. 10 together with its paramagnetic and 
antiferromagnetic Brillouin zone boundaries. The MBZ 
has full tetragonal symmetry and, in a transverse field, 
all the states at its boundary retain Kramers degeneracy. 
The leading terms in the expansion of g± (p) are linear 
near the MBZ faces, quadratic near the edges, and cubic 
near the vertices. 



CePd 2 Si 2 and CeRh 2 Si 2 

The heavy fermion metal CePd2Si2 has a body- 
centered tetragonal structure of ThCr 2 Si2 type, shown 
in Fig. 11(a). It is isostructural to CeCu2Si2 - the 
first discovered heavy fermion superconductor [40] 
and CeCu2Ge2, an incommensurate antiferromagnet [41], 
that becomes superconducting above 70 kbar in a pres- 
sure cell [42]. 

Below about 10 K, CePd2Si2 orders antiferromagnet- 
ically as shown in Fig. 11(a), with wave vector Q = 
[55O] , and a low-temperature ordered moment of about 



0.7 fiB per Ce atom. Its Sommerfeld coefficient is en- 
hanced to about 100 mJ/K 2 -mol. Samples of the present 
generation show residual resistivity in the p£l-cm range 
[22]. Under hydrostatic pressure of 26 kbar, the Neel tem- 
perature drops to under 1 K and, in a pressure window 
of ±5 kbar around this value, superconductivity appears, 
with a maximum transition temperature of about 0.4 K 
[43] . Curiously enough, normal state resistivity near this 
pressure follows a temperature dependence, that docs 
not fit the p{T) = po + AT 2 temperature dependence 
of the Landau Fermi liquid theory, but instead behaves 
as p{T) ~ T 1 - 2 over more than a decade in temperature, 
between about 1 and 40 K [43]. 




FIG. 11: (color online). Geometry of CePd 2 Si 2 in real and in 
reciprocal space, (a) Tetragonal unit cell of CePd 2 Si 2 , show- 
ing Ce atoms and the orientation of their magnetic moments. 
A full sketch, showing Pd and Si atoms, is given in [45]. (b) 
The Brillouin zone boundary of paramagnetic CePd 2 Si 2 is 
shown by dashed lines. The shaded hexagonal prism is its 
antiferromagnetic counterpart. 

The unit cell of CePd2Si2, and its paramagnetic and 
antiferromagnetic Brillouin zone boundaries are shown 
in Fig. 11. By symmetry, the degeneracy manifold in a 
transverse field includes the two hexagonal faces of the 
MBZ boundary, one of which is shown by darker shading 
in Fig. 11(b), and the four segments, two of which are 
shown in black. Along these segments, which are a three- 
dimensional analogue of point S in Fig. 6(a), another 
sheet of the degeneracy surface crosses the side faces of 
the MBZ. According to de Haas-van Alphen (dHvA) ex- 
periments [44], several Fermi surface sheets cross the de- 
generacy surface. The leading term in the expansion of 
g±(p) around the hexagonal MBZ faces is linear. 

CeRli2Si2 is an isostructural relative of CePd 2 Si2, with 
a modestly enhanced Sommerfeld coefficient of about 23 
mJ/K 2 -mol. Between T N1 w 36K and T N2 ~ 25K, it 
develops Neel order with Q = [55O] [45, 46]. Magnetic 
structure below T/V2 has not yet been established unam- 
biguously [45, 46]. Both T)vi and T/V2 drop under pres- 
sure [47] and, in an extended pressure window above 5 
kbar, CeRh2Si 2 becomes superconducting at a T c with a 
maximum of about 0.5K [48]. Antiferromagnetic struc- 
ture of CeRh 2 Si2 between T^i and T/v2 coincides with 
that of CePd2Si2, as does the degeneracy surface in Fig. 
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11(b). According to [49], at least one sheet of the Fermi had residual resistivities well below /if2-cm, and residual 
surface of CeRfi2Si2 crosses the degeneracy surface or resistivity ratios of over a 100 [54, 57]. 
comes close to it. 



Other materials of interest 



Neodymium hexaboride 

Rare earth hexaborides RBg are an interesting fami- 
ly, whose members show diverse electron and magnetic 
properties. Of the (relatively) simple ones, LaBg is a 
diamagnetic metal, and SmBg is a mixed valence semi- 
conductor. Of the ordered materials, EuBg is a ferro- 
magnetic semi-metal, and CcBg is a heavy fermion metal 
with at least two ordered phases, whose nature remains 
to be elucidated after nearly forty years of research. 

Three members of the family: NdBg, GdBg, and PrBg, 
are antiferromagnetic at low temperature. In PrBg [50] 
and in GdBg [51, 52] alike, two different low-temperature 
antiferromagnetic states have been found. 
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FIG. 12: (color online). Neodymium hexaboride. (a) Crys- 
talline and magnetic structure of NdE>6: a CsCl structure with 
B octahedra replacing CI atoms, and Nd in place of Cs. The 
arrows show magnetic moments of neodymium atoms, (b) 
Dashed lines show the cubic Brillouin zone of paramagnetic 
NdBg. The shaded square prism inside it is the tetragonal 
Brillouin zone boundary in the antiferromagnetic state. Its 
darker face denotes the degeneracy plane. 



Neodymium hexaboride NdBg presents a simpler pic- 
ture: below about 8 K, it is a collinear type-I antiferro- 
magnet with ordering vector Q = [00^] and an ordered 
moment of about 1.74 /j,b [53]; antiferromagnetism dou- 
bles its cubic unit cell in the [0 1] direction, as shown 
in Fig. 12(a). Thus the cubic magnetic Brillouin zone 
reduces by half in the [0 1] direction, while keeping its 
other two dimensions intact, as shown in Fig. 12(b). 

In a transverse field, the Kramers degeneracy is pro- 
tected at the two faces of the MBZ boundary, one of 
which is shown by darker shading in Fig. 12(b). Accord- 
ing to de Haas-van Alphen measurements [54, 55] and to 
calculations [56] , at least one sheet of the Fermi surface 
crosses the degeneracy surface. Recently studied samples 



This subsection contains a brief discussion of other an- 
tiferromagnets, where symmetry may protect the degen- 
eracy of special electron states against transverse mag- 
netic field, giving rise to Zeeman spin-orbit coupling (1). 

Cuprate superconductors: Electron-doped 
cuprates such as Nd2- a; Ce :! ;Cu04±5 develop commensu- 
rate antiferromagnetic order in a wide range of doping 
[58], albeit with a modest staggered moment [59]. 
For such materials, Fig. 6(b) describes the paramag- 
netic and antiferromagnetic Brillouin zone boundaries. 
Angle-resolved photoemission experiments [60] on 
Nd2- a; Ce a ;Cu04±5 have found carriers in a vicinity of 
the MBZ boundary. In a transverse magnetic field, these 
carriers are subject to Zeeman spin-orbit coupling (1), 
provided antiferromagnetism in the sample is developed 
well enough. 

Recent observation [61] of magnetic oscillations in 
YBa 2 Cu 3 06.5 testifies to great progress in sample quality 
of cuprates. And the fact that this and other underdoped 
cuprates are, at the very least, close to commensurate an- 
tiferromagnetism, makes them an interesting opportunity 
to examine the effects of Zeeman spin-orbit coupling. 

Borocarbides RT 2 B 2 C with R = Sc, Y, La, Th, Dy, 
Ho, Er, Tm or Lu and T = Ni, Ru, Pd or Pt have been a 
subject of active research, driven by interest in interplay 
between antiferromagnetism and superconductivity [4]. 
At low temperatures, commensurate antiferromagnetism 
develops in a number of borocarbides (for instance, in 
i?Ni 2 B2C with R = Pr, Dy or Ho), often with a large 
staggered moment (« 8.5 fiB for Dy and Ho) [4]. Zeeman 
spin-orbit coupling (1) is present whenever a sheet of the 
Fermi surface crosses the degeneracy manifold, and suc- 
cessful growth of high-quality single crystals [62] makes 
these materials an interesting case to study. 

Organic conductors are an immense and ever grow- 
ing class of quasi-low-dimensional materials, that show 
virtually all known types of electron states, found in con- 
densed matter physics [8]. Antiferromagnetism appears 
in several families of organic conductors, and manifesta- 
tions of Zeeman spin-orbit coupling (1) are likely to be 
found in some of them. 

Unfortunately, so far nearly all of the information on 
magnetic structure of organic antiferromagnets has been 
coming from indirect probes such as magnetic susceptibi- 
lity measurements [63, 64] and resonance spectroscopies 
[65-67]. Neutron diffraction studies are hampered by a 
typically small ordered moment, and by the difficulties 
of growing large enough single-crystalline samples. At 
the moment of writing, I am aware of only a single cycle 
of neutron scattering experiments [68-70] on an organic 
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conductor. Moreover, in families such as (TMTSF)2 
Bechgaard salts [8] and k-(BEDT-TTF) 2 X salts [71, 72], 
antiferromagnetic states are insulating, and their con- 
trolled doping remains a challenge [73]. 

With this word of caution, a number of organic conduc- 
tors may deserve attention. Semi-metallic Bechgaard salt 
(TMTSF) 2 N0 3 [74], developing a spin density wave state 
below about 9 K, may be one interesting case. Recently 
synthesized ethylenedioxytetrathiafulvalenoquinone-1,3- 
diselenolemethide (EDO-TTFVODS), that appears to 
turn antiferromagnetic below about 4.5 K, and remains 
normal down to the lowest studied temperature of 0.45 
K [75] , may be another. Finally, recent studies [76-78] of 
[Au(tmdt) 2 ], where tmdt denotes trimethylenetetrathi- 
afulvalenedithiolate, draw attention to this organic con- 
ductor. Albeit the material is not yet fully characterized, 
and its large single crystals remain difficult to grow, it ap- 
pears to have a Neel temperature of about 110 K [77, 78], 
which is anomalously high for an organic material - and 
shows normal conduction down to at least 10 K. 

Heavy fermion materials: Several heavy fermion 
antiferromagnets were reviewed in detail above. A num- 
ber of other interesting examples may be found in [10]. 

Gadolinium antiferromagnets (sec [79, 80], and 
Table I in [81]) offer two important advantages for an 
experimental study of the Zeeman spin-orbit coupling. 
Firstly, their often elevated Neel temperature Tjv (such 
as 134 K for GdAg, or 150 K for GdCu) facilitates ex- 
perimental access to temperatures well below Tjv, where 
thermal fluctuations of antiferromagnetic order are frozen 
out. Secondly, large ordered moment of these materials 
(about 7.5/ib for GdAg, and about 7.2/j.b for GdCu 2 Si 2 ) 
quenches quantum fluctuations. Therefore, gadolinium 
antiferromagnets fit well into the present framework with 
its neglect of both quantum and classical fluctuations - 
and shall be convenient for a study of various effects of 
the Zeeman spin-orbit coupling. 

Iron pnictides have been attracting immense atten- 
tion [82] due to appearance of commensurate antifcrro- 
magnetism [7] and high-temperature superconductivity 
[83] in this copper- free family of materials. Combination 
of commensurate antiferromagnetism [7] with essentially 
metallic normal state conduction [83, 84] not only con- 
trasts iron pnictides with the cuprates (that are believed 
to be Mott insulators), but also makes the former mate- 
rials likely to manifest a substantial momentum depen- 
dence of the (/-tensor. 

V. DISCUSSION 
Effects of relativistic spin-orbit coupling 

The arguments above appealed to the exchange sym- 
metry approximation [17]: the point symmetry opera- 
tions of the electron Hamiltonian in an antiferromagnet 



were considered inert with respect to spin, and the com- 
mon relativistic spin-orbit coupling, that appears in the 
absence of an external magnetic field, was thus neglected. 
I will now examine the effects it may have. 

Firstly, this spin-orbit interaction generates magnetic 
anisotropy, that creates a preferential orientation of the 
staggered magnetization n with respect to the crystal 
axes. In an experiment, this allows one to vary the mag- 
netic field orientation with respect to n as long as the 
field remains below the reorientation threshold. 

At the same time, the spin-orbit coupling may elimi- 
nate those spatial symmetries, that rotate the magneti- 
zation density with respect to the lattice. For instance, 
certain spin rotations and spatial transformations, that 
were independent symmetries within the exchange sym- 
metry approximation, may survive only when combined. 
I will now illustrate this by two examples of Section IV. 

A simple case of the spin-orbit coupling affecting the 
Kramers degeneracy manifold in a transverse field is 
given by a two-dimensional antiferromagnet on a square- 
symmetry lattice as in Fig. 6(b). Here, the Kramers de- 
generacy at the antiferromagnetic Brillouin zone bound- 
ary relies, everywhere except for points E, on the sym- 
metry with respect to reflections in diagonal planes 1 and 
2. If either of these reflections changes the orientation of 
A r with respect to the crystal axes, a spin-orbit coupling 
may lift the degeneracy at a relevant part of the magnetic 
Brillouin zone boundary, except for points S. However, 
if the magnetization density points along one of these di- 
agonal axes, the degeneracy survives at the two faces of 
the MBZ boundary, that are normal to this axis. 

In the case of commensurate order in chromium, con- 
sider a single-domain sample with magnetic structure 
shown in Fig. 7. For a Bloch state |p) with a momen- 
tum p at one of the two horizontal faces of the magnetic 
Brillouin zone in Fig 8(a), the degenerate partner state 
0T a U n (Tr)\p) has momentum — p at the other horizon- 
tal face of the MBZ. Coordinate rotation by n around 
the vertical symmetry axis, passing through the center T 
of the Brillouin zone, transforms the momentum p into 
p + Q, equivalent to p up to reciprocal lattice vector Q 
of the antiferromagnetic state. 

By contrast, for a momentum p at one of the vertical 
faces of the MBZ, coordinate rotation by n around a 
horizontal axis is required, and such a rotation inverts 
A r once the latter is attached to the crystal axes. Thus 
spin-orbit coupling tends to lift the degeneracy at the 
vertical faces of the MBZ, leaving it intact at the two 
horizontal faces. The other examples of Section III can 
be analyzed similarly. 

Finally, those spin-orbit coupling terms, that act di- 
rectly on the electron spin and tend to lift the double 
degeneracy of Bloch eigenstates even in the absence of 
magnetic field, were neglected here altogether. 
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Relation to earlier work 

When symmetries of a system involve time reversal - 
alone or in combination with other operations - a proper 
treatment must involve non-unitary symmetry groups: 
those containing unitary as well as anti-unitary elements. 
In this case, construction of irreducible representations 
is complicated by the fact that anti-unitary elements in- 
volve complex conjugation. In a group representation, 
combination of two unitary elements u\ and 112 is rep- 
resented by the product of the corresponding matrices 
D(«i) andD(u 2 ): D(wiu 2 ) = D(ui)D(u 2 )- By contrast, 
combination of an anti-unitary element a with a uni- 
tary clement u involves complex conjugation: D(au) = 
D(a)D*(u). As a result, irreducible representations of a 
non-unitary group must include a unitary representation 
and its complex conjugate on an equal footing. Discus- 
sion of such representations (called co-representations) 
was given by Wigner, along with the analysis of aris- 
ing possibilities with the help of the Frobenius-Schur cri- 
terion [85]. Later, Herring studied spectral degenera- 
cies, emerging in crystals due to time reversal symme- 
try and, among other things, extended this criterion to 
space groups [86]. In a subsequent work, Dimmock and 
Wheeler generalized the criterion further, to magnetic 
crystals, and pointed out the sufficient condition (6) for 
the appearance of extra degeneracies [13]. 

The present work identifies the symmetry, that pro- 
tects the Kramers degeneracy in a Neel antiferromag- 
net against transverse magnetic field, as a conspiracy be- 
tween the anti-unitary symmetry U n (7r)T a #, inherent to 
any collinear commensurate antiferromagnet in a trans- 
verse field, and the crystal symmetry of those special mo- 
menta at the MBZ boundary, that are defined by Eqn. 
(6). Formally, the present work is an extension of [13], 
since one may think of the last two terms in (4) as of the 
exchange field of a fictitious magnetic crystal in zero field. 
However, Kramers degeneracy in a magnetic field has 
rather special and remarkable experimental signatures, 
some of which are outlined at the end of this section. 

Last but not the least, Ref. [14] was an important 
source of inspiration for the present work. Its authors 
studied the electron eigenstates in a Neel antiferromag- 
net on a lattice of square symmetry and, for this particu- 
lar case, pointed out the disappearance of g± (p) at the 
MBZ boundary, as well as the ensuing substantial mo- 
mentum dependence of g± in the Zeeman coupling (1). 
The present article builds on Ref. [14] by elucidating the 
structure of the manifold of degenerate states for an ar- 
bitrary crystal symmetry, and for an arbitrary transverse 
field that can be sustained by the antiferromagnet before 
its sublattices collapse. This is to be contrasted with the 
analysis of Ref. [14], performed to the linear order in the 
field. Several other aspects of Ref. [14] are discussed in 
Appendix D. 



Experimental signatures 

The Kramers degeneracy at special momenta on the 
MBZ boundary and the resultant Zeeman spin-orbit cou- 
pling have a number of interesting consequences. For in- 
stance, a substantial momentum dependence of g±(p) in 
Eqn. (1) means that, generally, the Electron Spin Re- 
sonance (ESR) frequency of a carrier in the vicinity of 
the degeneracy manifold varies along the quasiclassical 
trajectory in momentum space. 

For a weakly-doped antiferromagnetic insulator with 
a conduction band minimum on the degeneracy mani- 
fold, this leads to an inherent broadening of the ESR line 
with doping and, eventually, complete loss of the ESR 
signal. In fact, this may well be the reason behind the 
long-known 'ESR silence' [87] of the cuprates. Suppres- 
sion of Pauli paramagnetism in the transverse direction 
with respect to staggered magnetization is another sim- 
ple consequence of vanishing g±(p). 

At the same time, a momentum dependence of g±(p) 
allows excitation of spin resonance transitions by AC 
electric rather than magnetic field [88, 89] - a vivid ef- 
fect of great promise for controlled spin manipulation, 
currently much sought after in spin electronics. Its ab- 
sorption matrix elements are defined by H p of Eqn. (19). 
Comparison with Eqn. (11) shows that, within the weak- 
coupling model (8), S p /ft is of the order of the antiferro- 
magnetic coherence length £ ~ ' an( ^ ma y be of the 
order of the lattice period or much greater. By contrast, 
the ESR matrix elements are defined by the Compton 
length Ac = ~ 0.4 pm. Thus, matrix elements of 
electrically excited spin transitions exceed those of ESR 
by about |§ • ^ w 137 • or at least by two orders 
of magnitude. Being proportional to the square of the 
appropriate transition matrix element, resonance absorp- 
tion due to electric excitation of spin transitions exceeds 
that of ESR at least by four orders of magnitude. 

Last but not the least - according to Eqn. (19), reso- 
nance absorption in this phenomenon shows a non-trivial 
dependence on the orientation of the AC electric field 
with respect to the crystal axes, and on the orientation 
of the DC magnetic field with respect to the staggered 
magnetization. 

The Zeeman spin-orbit coupling may also manifest it- 
self in other experiments on antiferromagnetic conduc- 
tors. In particular, de Haas- van Alphcn oscillations [90] 
and magneto-optical response may be modified. In va- 
rious types of electron response, interesting effects may 
arise due to an extra term ~vzso i n the electron veloci- 
ty operator, emerging due to a substantial momentum 
dependence of gj_(p) in Eqn. (1): 

vzso = V p Hzso = -AisV p .g ± (p)(H ± • a). (20) 

This term describes spin current. However, ,gj_(p) is even 
in p due to inversion symmetry and thus, in equilibrium, 
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the net spin current must vanish. This may change, if the 
system were tilted, say, by electric current or otherwise 
- however, the resulting effect would be proportional to 
the 'tilt' and, in addition to this, would be small in the 
measure of H±/A. 



Conclusions 

In this work, I studied the degeneracy of electron Bloch 
states in a Neel antiferromagnet, subject to a transverse 
magnetic field - and described the special points in mo- 
mentum space, where the degeneracy is protected by a 
hidden anti-unitary symmetry. 

I discussed the simplest properties and some of the 
manifestations of the Zeeman spin-orbit coupling, arising 
in a magnetic field due to this degeneracy, and outlined 
several examples of interesting materials, where such a 
coupling may be present. Finally, I reviewed the results 
and their relation to earlier work. 

The degeneracy of special Bloch states in a transverse 
field hinges only on the symmetry of the antifcrromag- 
netic state, and thus holds in weakly coupled and strongly 
correlated materials alike - provided long-range antifer- 
romagnetic order and well-defined electron quasiparticles 
are present. Under these conditions, thermal and quan- 
tum fluctuations of the antiferromagnetic order primarily 
rcnormalize the sublattice magnetization, leaving intact 
the degeneracy of special electron states in a transverse 
field - certainly in the leading order in fluctuations. De- 
tailed account of fluctuations is outside the scope of this 
article. 
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APPENDIX A: ORTHOGONALITY RELATION 



This Appendix proves the relation 

(0| [O0] + | [OO] ty) = <V#>, 



(21) 



where \<j>) and \ip) are arbitrary states, O is an arbitrary 
unitary operator, and is time reversal. In the main 



text, this relation is used for \<j>) = 09\ip); in this case, 
when read right to left, Eqn. (21) yields 

(VI00IV) - (H(O0) + ] 2 \(O8)\iP). (22) 

Whenever is an eigenvector of the linear operator 
[09] 2 with an eigenvalue different from unity, Eqn. (22) 
proves orthogonality of \ip) and 08\ip). 

The proof of Eqn. (21) is based on the obvious relation 
(C(j>,Cip) = (ip, </>) for arbitrary complex vectors <j> and tp, 
where (ip,4>) = X^V'T'^j denotes scalar product, and C 
is complex conjugation. Hence, for an arbitrary unitary 
operator O, one finds (OC(f>, OCip) — (ip,(t>), due to in- 
variance of scalar product under unitary transformation. 
Time reversal 9 can be presented as a product of C and 
a unitary operator [85]: 9 = VC , thus C = V~ 1 9 and, 
therefore, (O0<p, 09%jj) = (ijj, </>). As a result, for arbitrary 
states \ip) and \<j>), one finds (<j>\ [09} + | [OB] \ip) = (ip\<j>), 
which indeed amounts to (21). 



APPENDIX B: CANTING OF THE 
SUBLATTICES 

Canting of the two sublattices by transverse field Hj^ 
induces a component A^r of the magnetization density 
along the field, with the periodicity of the underlying 
lattice: A^ +& (H±) = A^(H±), as shown in Fig. 3. As 
a result, the diagonal part of Hamiltonian (8) acquires 
an additional term (Ap- • er), and Hamiltonian (8) thus 
takes the form 



n 



6 P - (Ap- • a) (A,, • a) 

(Ayo-) e p+Q - (Ap- +Q ■ a) 



(23) 



where Ap- 



H± + Ap-. 



The same choice of spin axes as in Section III splits 
Hamiltonian (23) into two independent pieces 



Wi( 2 ) 



An 



e p+Q 



± A 



p+Q 



(24) 



As in Section III, momentum shift by Q maps TL\ and TL2 
onto each other, and spectral symmetries of Hamiltonian 
(23) coincide with those discussed in the second subsec- 
tion of Section III. Thus all of the conclusions of Section 
III remain valid after sublattice canting is accounted for. 

However, while degeneracy at special points is pro- 
tected by symmetry, the shape of the manifold of de- 
generate states may change under various perturbations. 
For instance, sublattice canting in a transverse field mo- 
difies the equation, describing this manifold and, for the 
conduction band, turns it into 



Cp 



^PSp 



(25) 
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where P = §[A£ + A^ +Q ] and ^ p = §[A£ - A£ +Q ]. 
Since has the real-space periodicity of the paramag- 
netic state, Ap- enjoys the same reciprocal space symme- 
try as e p . In particular, Ap +2 q = Ap-, and Ap- = A^ p 
(the latter property is also protected by the U\8 symme- 



try). At the same time, ip p +Q 



-tp p , and (j) p+Q 



thus the symmetry-dictated degeneracy points such as 
p = y explicitly belong to the manifold of Eqn. (25), as 
they should. 

In the limit of vanishing H±_ , Ap- is linear in the field: 

Ap = H±[l + Xp], where xp~ describes microscopic trans- 
verse susceptibility of the antifcrromagnct. Now one may 
expand Eqn. (25) to linear order in the field to obtain 
the following equation for the degeneracy manifold: 



Xp 



xp^Cp 



= 0, 



(26) 







where Xp = X P ±Xp+Q- Compared with equation ( p 
of Section III, the sublattice canting affects the degene- 
racy manifold already in zeroeth order in H±. 



APPENDIX C: DIMENSIONALITY OF THE 
DEGENERACY MANIFOLD 

The dimensionality of the degeneracy manifold in a 
transverse field is one less than that of the momentum 
space for simple reasons, that rely only on the symmetry 
of the antiferromagnctic state. According to Eqn. (3), 
zero-field Bloch eigenstates |1) = |p) and |2) = XT a 9\p) 
form a Kramers doublet at momentum p. Its splitting 
S£(p) in a transverse field Hj^ is given by 



8£(p) = 2^|y 12 (p)| 2 + - [Vn(p) - V 22 {p)]\ (27) 

where Vij{p) = («|(Hj_ • <r)\j). Magnetic field being uni- 
form, (Hj_ • cr) commutes with TT a ; it also changes sign 
under time reversal. Thus, V 22 (p) = — Vii(p)- At the 
same time, the off-diagonal matrix element Vi 2 (p) vani- 
shes identically: 

<p|(Hl • <r)2T a (9|p) = £>|(Hl • <T)|q)(q|ZT a %) = 
q 

= ^y 11 (p)<S pq (q|JT a 0|p) = y 11 (p)(p|JT a 0|p> = 0, (28) 
q 

where insertion of unity 1 = ^ q |q)(q| was used in the 
first line, uniformity of Hj^ in the second, and the final 
equality followed from Eqn. (3). Therefore, 



8£( P ) = 2\v n (p)\, 



(29) 



and, barring a special case, equation 8£ (p) = de- 
fines a (d — l)-dimensional surface of zero <?j_(p) in d- 
dimcnsional momentum space. The Kramers degeneracy 



subset contains, at the very least, the star of the momen- 
tum p = Q/2 [see Eqn. (6) and the subsequent discus- 
sion], and the (k • p) expansion [91] around these points 
shows, that they are not isolated, but rather belong to a 
(d — l)-dimensional manifold. The latter is continuous, 
with the obvious exception of d = 1 . 

Finally, notice that, according to (29), 8£ (p) is peri- 
odic with the antiferromagnctic ordering wave vector Q: 



S£(p + Q) = 8£(p), 



(30) 



thanks to Q being a reciprocal lattice vector in the anti- 
ferromagnetic state. Therefore, properties (16) and (17) 
are indeed model-independent, as opposed to hinging on 
an approximation of the weak-coupling model (8). 



APPENDIX D: REVISITING [14] 

In Rcf. [14], Brazovskii and Lukyanchuk stated that 
operator A = (n-er) exchanges the momenta p and p + Q 
in (8), and thus represents the momentum boost by the 
ordering wave vector Q in reciprocal space. In a com- 
mensurate antiferromagnetic state, Q becomes a recip- 
rocal lattice vector, and thus A must be a symmetry of 
the Hamiltonian. With the assumption of the effective 
Zeeman coupling (1), this lead the authors of Ref. [14] 
to the relation g±(p + Q) = — <7_i_(p) (Eqn. (16) of the 
present work), and to the conclusion that g±(p) = at 
the MBZ boundary. Unfortunately, while this beautiful 
result is indeed correct for a lattice of square symmetry, 
several circumstances prevent one from embracing these 
arguments. 

Most importantly, they hinge solely on the symme- 
try under translation by Q, put otherwise - on commcn- 
surability of magnetic order with the crystal lattice. If 
correct, this would imply that, in an arbitrary commen- 
surate Neel antifcrromagnct, Kramers degeneracy takes 
place at the entire MBZ boundary regardless of the un- 
derlying crystal symmetry. The toy example (18) with 
77 ^ 1 shows, that this is not at all necessarily the case. 

Indeed, for a generic crystal symmetry, the condi- 
tion g±(p + Q) = — <7_i_(p) (see [14] and Eqn. (16) of 
the present work) does not, by itself, restrict the ma- 
nifold <7_i_(p) = to the MBZ boundary. However, the 
'Kramers' subset of the manifold of degenerate states can 
be obtained by combining Eqn. (16) with the crystal 
symmetries (see the subsection on spectral symmetries in 
momentum space in Section III). For instance, combined 
with the inversion symmetry g±(— p) = <7_i_(p), Eqn. (16) 
stipulates that g±(Q/2) = 0. Similarly, disappearance of 
g±(p) at the entire MBZ boundary for the square sym- 
metry case can be obtained by using Eqn. (16) and the 
point symmetries of the square lattice. For a finite as op- 
posed to infinitesimal field, these results were established 
in Section II and in the first two examples in Section IV. 
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On a more technical level, the operator A is equiva- 
lent to U n (7r) and thus inverts the sign of the transverse 
component of the field. Hence, in a field with non-zero 
transverse component H^, A ceases to be a symmetry of 
the Hamiltonian, in agreement with the first line of Ta- 
ble I - and thus can no longer represent the momentum 
boost by Q. 
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